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Abstract 

The first in a long series of papers by John T. Lewis, in cohaboration with G. W. Ford and 
the present author, considered the problem of the most general coupling of a quantum particle 
to a linear passive heat bath, in the course of which they derived an exact formula for the free 
energy of an oscillator coupled to a heat bath in thermal equilibrium at temperature T. This 
formula, and its later extension to three dimensions to incorporate a magnetic field, has proved to 
be invaluable in analyzing problems in quantum thermodynamics. Here, we address the question 
raised in our title viz. Nernst's third law of thermodynamics. 
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I. INTRODUCTION 



In recent years, the realization that quantum effects can play an important role in thermo- 
dynamic theory has led to an intense interest in quantum and mesoscopic thermodynamics 

0; y ; 0; Q| ; to the cxtcut that some authors jl[ 3,0,0] have questioned whether the fun- 
damental laws of thermodynamics remain valid. The arguments are often very subtle and 
sometimes not easily dismissed. Here, our goal is not to survey the whole area but, instead, 
to focus in depth on questions which have been raised in relation to Nernst's third law of 
thermodyanmics, which states that the entropy S vanishes as the temperature T 0. Ford 
and the present author ^ already considered this question for the one-dimensional problem 
of a quantum oscillator in an arbitrary heat bath and also briefly commented on what we 
considered to be the flaw in the calculations of other authors. We now expand on these re- 
marks in Sec. II (where we discuss various approaches to the calculation of S), by carrying 
out a detailed calculation in which we demonstrate explicitly how use of approaches based 
on an incorrect utilization of the Wigner function formalism, and also the von Neumann 
formula, leads to an incorrect form of a well-known result. By contrast, the method we 

n 

employed in [5] was based on calculating S from an exact expression for the free energy F 
calculated by us in 1985 at the start of a long collaboration with John T. Lewis 0. This 



formula for F was later extended to three dimensions to incorporate a magnetic field 
Thus, in Sec. Ill, we review these results and explicitly write down the results for F and 
S in the three-dimensional case. All of these results involve just a single integral and, in 
Sec. IV, we evaluate this integral for the case of an Ohmic bath and low temperature. This 
enables us to confirm Nernst's law. Our conclusions are given in Sec. V. 

II. CALCULATION OF THE ENTROPY S 

For a quantum particle in a heat bath, a basic quantity is the density matrix 

p = e-^/^^/Tr{e-^/'=^}, (2.1) 

where H is the Hamiltonian for the whole system (quantum particle plus heat bath plus 
interaction). The question is how to calculate the entropy S from this expression. In Ref. 4, 
three different results for the entropy are obtained based on use of the free energy {Sp given 
in (3.59) of ^), the Wigner distribution {Sb given in (4.5) and (4.38) of Q| and referred to 



as the " Boltzmann" energy) and the von Neumann expression {S^n given in (4.36) of P|). 

In particular, the authors of ^ argue that " - - neither the von Neumann entropy nor the 

Boltzmann entropy vanishes when the bath temperature is zero" , leading them to note " - - 

the violation of the third law reported here for nonweak coupling" . On the contrary, we argue 

that all such calculations should lead to the same result so we are motivated to examine the 

question of what went wrong. In the process, we will argue that a calculation based on the 

free energy is the simplest way to proceed to obtain the correct result, particularly if one 

uses the formula given in Qj. 

First, we examine the calculation of S^n, which was based on applying the von Neumann 

formula to the reduced density matrix. But as we pointed out in p|, for a system with a 

non-negligible interaction energy, " - - the von Neumann formula can only be applied to 

the entire system and not to the reduced system" . At first glance, these may seem to be at 

variance with the fact that numerous calculations rely on the calculation of the von Neumann 

entropy of mixed states obtained by tracing out other variables. The answer is that this 

procedure is correct only when the interaction energy is negligible, as is often the case (so 

that the total entropy is equal to the sum of the entropies of the parts). However, this is not 

so in the case of heat baths (where we are dealing with non-additivity of entropies), as may 

be verified explicitly. This is what motivated our calculation with John Lewis, where the 

free energy F of the electron in the presence of the heat bath was calculated by first using 

the von Neumann formula to calculate the free energy of the entire interacting system and 

then subtracting the well-known free energy of the heat bath itself. Simply differentiating 

- F with respect to T [see [2.10) below] leads to the correct result Sp. The authors of |^ 

used our procedure in their calculation of Sp [their Eq. (3.59)], obtaining a correct result 
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in agreement with our previous calculation of this quantity [5], but in disagreement with 
their result [given in their Eq. (4.36)] for S^n (which was obtained assuming additivity of 
entropies). If, instead, one uses the density matrix of the whole system, the corresponding 
result will coincide with the result Sp. 

Second, we examine the calculation of Sb of Q] which was based on writing down the von 
Neumann formula for the entropy, S = — A;Tr{p logp}, and then making use of the reduced 
Wigner distribution function W |S] corresponding to p, by simply replacing p with W in the 
von Neumann expression, as may be verified by examination of equations (4.3) and (4.4) of 
0], leading to the authors' equation (4.5). The calculation within this framework is correct 
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but, unfortunately, the framework is not correct, as we shall now demonstrate explicitly in 
the simplest case where the interaction energy is taken to be negligible (i.e. by letting the 
dissipation parameter 7 — 0), where well-known results for the entropy of a single oscillator 
in equilibrium at temperature T [9] may be used as a benchmark. For an oscillator at 



temperature T, we have the well-known exact result for the Wigner distribution 

W{q,p) = Nexp{-aH), (2.2) 

where 

i^^coth(M), (,3) 

N=^ = (nh)-' tanh (^^) = (27r^)-^ N^, (2.4) 

^(g,p) = |^ + ^mwV- (2.5) 
J J dqdp W{q,p) = 1, (2.6) 

< Hiq,p) >= 1 1 dqdp Wiq,p) H{q,p) = ^. (2.7) 

If we now replace (incorrectly as we shall see) p in the von Neumann formula by 27r^ W 
then we obtain for the entropy (which we denote here as S^^\ indicating that it was derived 
using the Wigner distribution) 

{S^^^k) = - J J dqdpW{q,p)\og[27rhW{q,p)] 

= —\ogNi J J dqdp W{q,p) + a J J dqdp H(q,p)W{q,p) 
= —logNi—aN-^ f [ dqdpexp{—aH) 



and 



Also, we note that 



and 



da 

1 fi/jj 

log 



1 , hujp' 

- coth 



2 2 



+ 1. 

(2.8) 



This result is incorrect, being in striking contrast with the correct result given by 



(S/k) = - log [2 sinh(^c<;/2A;T)] + ^ coth{huj/2kT) . (2.9) 
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Only in the high temperature classical regime do equations (2.8) and (2.9) agree, in which 
case the Wigner distribution function simply reduces to a classical distribution function. 
One of the problems with this method stems from the fact that the Wigner distribution 
corresponding to log p is not the same as the log of the Wigner distribution corresponding 
to p. Thus, we need an alternative approach. 

In fact, as already noted , the best approach is simply to calculate F from the 

formula 



S{T) = - 




(2.10) 



and in the next section we discuss the method by which F is obtained. As we shall see, 
the .esu>t of ,4] . eo„eet a,.ee. w.th a .esu. p.ev.ous.. obtained us Q i„ tbe 

case of an Ohmic heat bath in the absence of a magnetic field. However, we will also note 
that the method used in for the calculation of F, is unnecessarily complicated and only 
applies to the Ohmic model whereas, by contrast, our result is expressed as a simple integral 
in terms of the generalized susceptibility and applies, not only for the Ohmic bath, but for 
an arbitrary heat bath and in the presence of a magnetic field. 



III. FUNDAMENTALS 

The Hamiltonian for a charged quantum oscillator (with a force constant K = mcug) 
moving in an external magnetic field and linearly coupled to a passive heat bath (consisting 
of an infinite number of oscillators) may be written as 

Ho = —\p--aY + -Kf^ + + -m,LU^(q,-fY]. (3.1) 

2mL^ c J 2 V^Sm,- ^ 2 ^ ^^^^ /J ^ ^ 

This is the independent oscillator (10) model in the presence of an external magnetic field 
B, where e, m, p, r are the charge, mass, momentum and position of the oscillator respec- 
tively and the corresponding quantities with the lower indices j refer to the jth heat bath 
oscillator. The vector potential A is related to the magnetic field B through the equation 

B{f) = V X A{r) . (3.2) 

Next, using the Heisenberg equations of motion, the problem has been formulated exactly 
in terms of the quantum Langevin equation. 

mf + f dt' p{t - t')'f{t') - -'r^B + Kf = F{t). (3.3) 

J —oo C 
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Here F(t) is the fluctuation force (wliose exact form is known but is not relevant to our 



present discussion) and /i(t) is the so-called memory (non-Markovian) term given by [7[ 



n{t) = '^mjtu'^ cos{uj jt) 9 (t), 



(3.4) 



where 6{t) is the Heaviside function. 

We note that fi{t) does not depend on the magnetic field. Fourier transforming (3.3), we 
obtain [3| 



where 

is the generalized susceptibility, with 

deW{u) 



A^^po- — \^~) BpBcr — tpariBrfXiuo- / detD{uo) 



A 



A2 _ (^1)2 

c 



and 



Also 



-1 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



/•oo 

= / dte'^^'fiit), (3.9) 

Jo 

/oo 
dte'^^'r^it), (3.10) 
-oo 

and so on, and where ep^ri is the Levi-Civita symbol, a totally antisymmetric tensor. 
The Greek indices stand for three spatial directions (i.e. p, a, etc. = 1,2,3) and we adopt 
Einstein summation convention for repeated Greek indices. From (3.4) and (3.9), we obtain 
the Fourier transform of the memory function 

1 1 



" \ ^ 2 



+ 



LU — UJ4 



(3.11) 



We have now all the tools required in order to calculate various quantities. Here, we discuss 
the free energy and associated entropy. 

The free energy ascribed to the oscillator, F{T), is given by the free energy of the system 
minus the free energy of the heat bath in the absence of the oscillator. This is a non-trivial 
quantity to calculate, details of which may be found in [3] leading to the result 

Fo{T,B) = - / dujf{uj,T)Im{—ln\deta{uj + io+) }, (3.12) 

TT Jo ^ duj L J 



where f{uj,T) is the free energy of a single oscillator of frequency u, given by 



f{uj,T) = kT\og[l-exp{-hu/kT)]. (3.13) 
Here the zero-point contribution {huj/2) has been omitted and 



det«(cu) = [a^^\uj)]^ 

with 



1 



eBio 



2, 



-1 



(3.14) 



[a(o)(cj)]-i = -mu^ + K-iuofi{uj). (3.15) 

This enabled us to write 

F,{T,B) = F,{T,0) + AF,{T,B), (3.16) 

where 

F(T,0) = - duj f(uj,T)Im{—lna'~^\uj)} (3.17) 
n Jo ^dui J 

is the free energy of the oscillator in the absence of the magnetic field (in agreement with 
Eq. (5) of Ref. 6, except for the extra factor of 3 which results from our consideration here 
of three dimensions) and the correction due to the magnetic field is given by 

AF,(T,i?) = -- rduf{uj,T)Im{^ln[l-{^f{a'^'\u)f]]. (3.18) 
TX Jo ^duo c 

IV. OHMIC HEAT BATH 

In the case of the Ohmic heat bath, fl{uj) = m% a constant, which is the simplest 
memory function one can choose and is an oft-used jj] choice. Thus, making use of (3.7) 
and (3.8), (3.3) becomes 

1 r°° 

F„(T,5) = - / dujf{uj,T) (4.1) 
7(^' + ^o) 7(^' + ^o) 7(^' + ^o) 1 



TX Jo 



\ (u;2 _ ^;2^2 _j_ ^2^2 (^^2 _ 1^2 _|_ a;^a;)2 _|_ ,^2^2 (^^^2 _ ^2 _ ij^ij^ -|- (^2^2 j ' 

where Uo = (K/m)^ and cUc = {eB/mc). 

We now wish to examine the low-temperature behavior of this result. First, we note 
that the function f{uj,T) vanishes exponentially for u ^ kT/h. Therefore as T — ^ the 
integrand is confined to low frequencies and we can obtain an explicit result by expanding 



the factor multiplying f{uj,T) in powers of u. Hence, in particular, we see that the terms 
involving the magnetic field (the Uc terms) are negligible since they always contain an lj 
factor. Thus our result corresponds to the results obtained in the absence of a magnetic 
field 0, except for the extra factor of 3 which results from our consideration here of three 
dimensions. Thus, proceeding as in j^], we obtain in the low-temperature case 

F{T) = r dujXog [1 - exp {-huj/kT)] 

tcluq Jo 

SO that 

Hence S{T) — > as T — 0, in conformity with the third law of thermodynamics. We also 
note that the result for the entropy is independent of the magnetic field B and actually 
corresponds (except for a factor of 3 because here we considered 3 dimensions) to the result 
obtained earlier by us 5] in the B = case and also agrees with the result for Sp given 
in 0]. The latter calculation was also based on a determination of the free energy F but 
this calculation was much more complicated, involving a non-trivial determination of the 
uhe interacting system By contrast, our calculation was based on a 



frequencies of t 
simple integral 



0, 



given here by (3.12), and which involves only the specification of the 



generalized susceptibility a{u). 
V. CONCLUSION 

Similar results to those obtained in the Ohmic case may be obtained also in the case of 
a blackbody radiation heat bath and, in fact, in the case of arbitrary heat baths. Thus, we 
conclude that for the case of a quantum oscillator coupled to an arbitrary heat bath, we 
have shown that Nernst's third law of thermodynamics is still valid: the entropy vanishes at 
zero temperature. In this connection we should emphasize that the basis of our discussion 
is the remarkably srmple formula for the free e,rerg,v of arr osc.Uator. in a,r arbitrary heat 
bath at arbitrary temperature T, obtained in the 1985 paper p| with J. T. Lewis. Since the 
validity of the second law of thermodynamics has also been called into question [l|, Q, 0, 0] , 
we note that this formula for the free energy was also used in refuting speculations that 



quantum effects could lead to extraction of energy from a zero-temperature heat bath 11 1. 
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